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Abstract 

Multitime evolution PDEs for Rayleigh waves are considered, using 
geometrical ingredients capable to build an ultra-parabolic-hyperbolic 
differential operator. Their soliton solutions are found based on ap- 
propriate hypotheses and Bernoulli ODEs. These multitime solitons 
develop complex behavior of deformation phenomena. 

Section 1 presents the single-time Rayleigh wave equations. Sec- 
tion 2 analyzes the geometric characteristics (fundamental tensor, lin- 
ear connection, vector fields, tensor fields) associated to the Rayleigh 
PDEs, showing the existence of an infinity of geometrical structures 
such that the multitime PDEs are prolongations of single-time PDE. 
Section 3 builds the ultra-parabolic-hyperbolic differential operator 
defining the multitime Rayleigh wave equations. Section 4 gives the 
technique capable to produce the multitime Rayleigh solitons. Sections 
5 and 6 praise explicit formulas for the multitime Rayleigh solitons. 

Mathematics Subject Classification 2010: 37K40, 35C08, 35Q51. 
Key words: multitime solitons, single-time Rayleigh PDE, multitime 
Rayleigh PDE, Bernoum ODE. 

1 Single-time Rayleigh PDEs 

In the physical and mathematical literature [l]-[4], [7]- [9], we find the Rayleigh 
wave equation 

utt - Uxx = e{ut - ul) (R) 
related to Rayleigh wave equation of Van der Pol type 

Utt - Uxx = e(l - u^)ut. (RVP) 
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Each of these has been used to model physical phenomena. Now, the PDE 
(R) serves as a model for the large amplitude vibrations of wind-blown, 
ice-laden power transmission lines, in time that, the PDE (RVP) describes 
plane electromagnetic waves propagating between two parallel planes in a 
region where the conductivity varies quadratically with the electric field. 

Just as their counterparts from ordinary differential equations, the PDEs 
(R) and (VP) can be transformed one to another. Their solutions can be 
obtained by simple operations performed on the solution of a certain first 
order, nonlinear wave equation. 

An initial-boundary value problem for Rayleigh nonlinear wave equation 
can be considered to be a simple model to describe the galloping oscilla- 
tions of overhead power transmission lines in a windfield. One end of the 
transmission line is assumed to be fixed, whereas the other end of the line 
is assumed to be attached to a dashpot system. 

Rayleigh surface waves are of particular importance in seismology, acous- 
tic, geophysics and electronics applications. 

Remark Some papers [6] describe as Rayleigh wave a type of seismic 
surface wave that moves with a rolling motion that consists of a combination 
of particle motion perpendicular and parallel to the main direction of wave 
propagation. The amplitude of this motion decreases with depth. Like pri- 
mary waves, Rayleigh waves are alternatingly compressional and extensional 
(they cause changes in the volume of the rocks they pass through) . 

2 Multitime PDE systems 

The passing from systems of PDEs with a single-time variable t to related 
PDE systems with m > 2 evolution variables t = {t""), a = l,...,m, is 
substantially complicated due to necessity of praising some reasons and some 
techniques of such change. The most natural way of changing is to use 
geometrical ingredients (derivation, trace etc) that extend the initial system. 
The theory and a systematic procedure for the construction of such new 
systems is presented here in the context of Rayleigh nonlinear waves. 

This paper provides new results regarding the multitime solitons in two 
and more temporal dimensions that can be of interest in physics. We over- 
pass the complexity and, furthermore, the difficulty of performing hand 
computations for Rayleigh PDE systems involving many temporal variables 
by using the symbolic software package in MAPLE. Our source of inspira- 
tion for introducing and studying multitime soliton PDEs is the paper [5]. 
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Also the papers [10]-[17] contains a lot of ideas in this direction. 



3 Multitime geometrical prolongations of 
Rayleigh PDEs 

Let us introduce and study some multitime geometrical prolongations of the 
Rayleigh PDEs, using related connection, fundamental tensor field, vector 
fields, tensor fields which leave on the jet bundle as ingredients in the Dif- 
ferential Geometry of the manifold RJ^. 

Suppose the multitime t = {t^, € IK*" is a parameter of evolution. 

We endow the manifold (jet bundle of order one) J^(1R™ x IR, IR x IR™) 

/ ()u\ 

with a distinguished symmetric linear connection F^g = F^g it,x,u,—j, 

and with a distinguished fundamental symmetric contravariant tensor field 
/ du \ 

h = I h'^^{t, X, u, — ) I of constant signature (r, z,s), r + z + s = m. Using 



dt , 

a function u : IR*" x IR ^ IR, we build the Hessian operator 

d'^u du 
{Hessru)afi = - Kpg^' «, /3, 7 G {1, 

and its trace, called ultra-parabolic-hyperbolic operator, 

We define a multitime PDE as 

□r,ft«-^ = 0, (1) 

where t = if} , f^) G i?"* and x ^ R. 
- Let 

C''{t,x,r],C), 7 = l,...,m 
be a distinguished vector field and 

be a distinguished tensor field. If we adopt the hypothesis 

h'^^it, X, rj, Orlsit, X, rj, 0^7 = V, 0^7 " ^"^^(*' O^a^0^i, (2) 
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then we obtain the multitime Rayleigh PDE 

.ad d^y^ du g.,. du du du d^u _ 

dt^dtP dt^ dt^dtPdt^ ^- w 

- If C^{t, X, r], 7 = 1, m and D'^{t, x, rj, 7 = 1, m are distin- 
guished vector fields and the constraint relation is 

h-^{t, X, ri, Orl/, it, X, ri, 0^7 = C^t, X, ri, ' D-'iP, V, OvH-y, (2') 
then we get a multitime Rayleigh wave equation of Van der Pol type 

""^^-'-l^-^-l^-SS-. (3') 

Theorem 

(1) There exists an infinity of geometrical structures T^^, h"^ , , B"'^'^ 
on such that a solution of the Rayleigh PDE (R) is also a solution of 
the multitime Rayleigh PDE (3). 

(2) There exists an infinity of geometrical structures T^^^, /i"'^, C"' , D'^ 
on such that a solution of the Rayleigh PDE (RVP) is also a solution 
of the multitime Rayleigh PDE of Van der Pol type (3'). 

Proof Let t^ = t and u = u{t^, x). 

(1) Suppose u = u{t^,x) is a solution of single-time Rayleigh PDE (R). 
The function v{t^, f^, x) = u{t^, x) is a solution of the multitime Rayleigh 
PDE (3) if the family of geometrical structures F^^^, /i"'^, C^, B"'^'^ is fixed 

by 

It is obvious that we have an infinity of geometrical structures that satisfy 
this algebraic equation. 

(2) Suppose u = u{t^,x) is a solution of single-time Rayleigh PDE 
(RVP). The function t™, x) = u{t^,x) is a solution of the multi- 
time Rayleigh PDE (3') if the family of geometrical structures F^^^, h"^ , 
CT, is fixed by 

h-^Tl^C, = C%-D\^^,. 

It is obvious that we have an infinity of geometrical structures that satisfy 
this algebraic equation. 

The foregoing Theorem justifies the term multitime geometrical prolon- 
gations of the Rayleigh PDEs. 
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4 Multitime Rayleigh solitons 



The first aim of this Section is to find some multitime solitons solutions 
for this multitime PDEs. In spite of the mathematical beauty, the distance 
between theoretical multitime models and real situations where they apply 
is far from our understanding. 

Let (/):/c]R— s-lRbea function of class C^. We seek for solutions of 
the PDEs (3) and (3') in the form of multitime solitons 

U{x,t) = ct>{x-Xat'') = ct>{z), (4) 

where Aq, a = 1, ...,m, is a constant vector and z = x — Xat". Then, the 
partial derivatives of the unknown function u{x, t) are 

S = *"W. ^ = ^ = (5) 

Substituting these derivatives in the PDEs (3) and (3'), we obtain the second 
order ODEs, 

[^"^A„A^ - 1] (p'\z) - 5"^TA„A^A^ (p'izf + C^A^ (j)'{z) = (6) 
and respectively 

[h'^^XaXff - 1] <f>"{z) - D^X^ (p\z)(l)'iz) + C^X^ (j)'{z) = 0. (6') 
Summarizing, we have 

Theorem If (f){z) is a solution of second order ODE (6) or (6') , then 
u{x,t) = (l){x—Xat°') is a multitime soliton solution of the multitime Rayleigh 
PDEs (3) or (3') respectively. 

In order to find some multitime Rayleigh solitons, we make some par- 
ticular choices of the elements that appear in the construction of the two 
multitime PDEs. 

First, for multitime Rayleigh solitons, wc consider the ODE (6) and we 
relate the metric tensor h'^^(t,x,r],^), the tensor field B'^^'^{t,x,r),£,), the 
vector field C^{t,x,r],^), 7 = l,...,m and the constant vector Aq, by the 
conditions 

' h^f'it, x, r], OXaXp -l = a{x- A«t°) = a{z) / 0, 
< X, V, OAaA^A^ = b{x - A„t°) = b{z), 

C'^{t, X, 77, ^)A^ = c{x — Aat") = c{z). 
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With these conditions, the ODE (6) becomes 

aiz)(l>"-biz)i(f>'f + c{z)(f>' = 0. (7) 

Second, for multitime Rayleigh sohtons, we take de ODE (6') and we 
relate the metric tensor h"l^{t,x,'r],^), the vector fields C'^{t,x,ri,(^) and 
D^(t,x,r],^), 7 = 1, ...,m and the constant vector Aq by the conditions 

' h^^{t, X, ri, 0AaA;3 - 1 = a(x - A^i") = a{z) 7^ 0, 

< D^{t, X, T], ^)A^ = b{x - A„r) = d{z), 
C"'{t,x,r],$,)X^ = c{x — Xat"') = c{z). 

With this choice, the ODE (6') becomes 

a{z)(t>" - d{z)(t)^<p' + c{z)(^' = 0. (7') 

We are looking for some solutions of the ODE (7) and (7'), with a view 
to finding multitime Rayleigh solitons. 

5 Families of multitime Rayleigh solitons 
5.1 Case of coefficients depending on z 

Denoting =* ■0, the second order ODE (7) becomes a first order ODE, 

a(z)V' - + c(z)V = 0, (8) 

called Bernoulli ODE. The general form of this ODE is called Abel ODE of 
the first kind, and it arose in the context of the studies of Niels Henrik Abel 
on the theory of elliptic functions, and represents a natural generalization 
of the Riccati equation. 

Since a{z) 7^ 0, we write the equation in the form 

a{z) a{z) 

By a change of the unknown function, ^ = tj^"'^, the Bernoulli ODE becomes 
a linear ODE, 

{'-244e = -2M, (10) 

a{z) a{z) 
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with the solutions 

c{z) 



a(z) 



dz 



K 



J aiz) 



ciz 



dz 



«(^) dz 



, K eM. (11) 



Since V' = C 2 , we obtain 

c{z) 



aiz) 



dz 



K -2 



[M. 

J aiz) 



{z) 



'-^dz ■ 
"(^) dz 



that is 



^\z) = 



c{z)_ 
a{z) 



dz 



a[z) 



-2 / -^dz 
«(^) dz 



(12) 



Therefore, we have found solutions of the multitime PDE Rayleigh (3): 
Theorem // we fix the above elements by the conditions 

[t, X, V, Aa A;3 - 1 = a{x - Aai«) 7^ 0, 

B"^^(t, x, rj, OAaA^A^ = b{x - A«t"), 
C^(t, X, rj, ^)X^ = c{x - Xat"), 



then the function 



u{x,t) = (t){x — Xat") 



represents a multitime soliton- solution for the multitime PDE Rayleigh (3) 
for every (p given by 



c{z)_ 
a{z) 



dz 



\ 



K-2 



-2 / --4dz 



dz, K eM. 



(13) 



Me V a{z)-^^ 
a{z) 
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5.2 Case of constant coefficients 

If we fix the elements h,B,C,X by the conditions 

h°'^{t, X, r], 0^Q^,3 ~ 1 = constant = a 7^ 0, 
< B°''^'^{t,x,ri,S^)XaX[sX'y = constant = b, 
C^{t, X, r],(^)X^ = constant = c, 

then the ODE (7) takes the form 

a0" - b{4>'f + c0' = 0. 

Denoting cp' ^= ip, the second order ODE with constant coefficients (14) 
becomes a first order ODE, called Bernoulli ODE 



(14) 



at/;' - bip^ + cV' = 0. 
The ODE (15) has separable variables, i.e., it can be written 

dip a f dip 



bipfi — cijj 



dz <^ 



I 



dz 



and we get the next sequence of equivalences 



dz 



<^ In 



-z + l, leJR 



(15) 



that is 



<^ ip^ - - = pip^ea'^, peJR* 
b 



^iP{z)=±. 



b{l-pe^') 



-, p e JR*, with 



— > 0, 

b{l-pe^') 



dz, p e M*, with > 0. (16) 

6(1 -pe~^) 6(1 -pe~^) 



8 



To calculate the primitive from the right hand side, we amplify the fraction 
by and then make the change of variables = t. We obtain the 
integral 



dt 



and a new change of variables, ~^ — gives a new integral. 



c J y h{s'^ — p) 



ds. 



There are two separate possibilities: 

c 

a) If we take ^ > 0) then 



J = ± 



a c 



b J s/s^ —p 



ds. 



- For p > 0, the integral is 



J = ± 



\/?ch-^ 


s 


V 





that is 





1 


V 





+ r, r G IR, 



+ r, reJR, peM* 



and we can write 



(p{z) = ±^ J '^ch-^ [Ke a^]+r, r G IR, K eM* 



c 
—z 



- For p < 0, we have 



J = ± 



It follows 



that is 



+ r, r G IR, 
+ r, r € IR, € IR+ 



J = ±-. \ sh-i 
c V 





1 







(17) 
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and then 

(z) = ±^^^ sh-^ {Ke~a^ \ +r, reJR, if G M;. (18) 

c c 
b) If we suppose 7 < 0, then, via 5^ — > 0, it foUows p > and 

b 6(1 -pe-^) 

the integral J becomes 

a f 1 

J = ±-\ hr / / o ds, p>o 

c V b J ^p- 

a I — c I s \ 

<^ J = ±-\ — arcsin — + r, r G IR, G Ml. 



c V 6 

It follows 



(f>{z) = ± - W — ^ arcsin ( ^ 
c V 6 yea , 



+ r, r G IR, p G Ml, 

VpJ 



that is 



0(z) = ± arcsin ^ivTe j + r, r G IR, if G IR+. (19) 

Therefore, we have found three families of solutions of the multitime 
Rayleigh PDE (3): 

Theorem a) If we fix the above elements by the conditions 
h'^^XaXg = constant ^ „ , — 7 — — > 0, 
then we get two families of soliton solutions 

( C-rx, ^ 



<x, t) = ± — W , , ch 



c^Xj y B'^f^^XaX/sX^ 



V 



(20) 
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where r G IR, KG ]R^, and 
u{x, t) = lb 



w/iere r G IR, ii' G 
b) if we take 



5°/5tA„A«A 



sh 



Ke 



h^n^Xfs - 1 



h'^^XaXa = constant, 



C^X^ 



B'^f^^XaX^X^ 
then we have another family of soliton solutions, 

h'^^XaXs-l I 



<0, 



u{x, t) = ziz 



-C^X^ 



X arcsin 





C^X^ 




C^X^ 


Ke 


h^^XaXp 


\ 





B'^P^XaXpX^ 



-(x-A«r)^ 



(21) 



(22) 



where reM, K G IR+. 

5.3 Mac-Laurin series soliton of 
multitime Rayleigh PDE 

In a previous section, we obtained the equation (7), 

a{z)(t>" -h{z){<t>'f + c{z)<t>' = 0, 

a second order ODE in (j). One assume that it has a solution which is 
analytic on an interval around z = and we search a Mac-Laurin series 
solution. Then we express as a power series in the form 



n=0 



(23) 



and we try to determine what the a^'s need to be. The resulting power 
series need to converge on an interval around origin. 
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We compute (t)'{z) and ^'^(2;) and (t)"{z): 

00 00 
ct>'{z) = J2 nan = E + = E A 



00 

n=l n=0 n=0 



00 / n 



(<^0'(^) = E E/3fc/3n-fc 
n=0 \jfc=0 / 



00 / n 



E E«fc+i "n-fc+i(A; + l)(n - A; + 1) = E 7n -z", 

n=0 \fe=0 / n=0 



00 / n \ 



(0'f(z) = E E7^/^n 
n=0 \i=0 



00 / n I 



= E I E E "'fc+i "n-j+i(A; + 1)(^ - A; + l)(n - z + 1) 1 z", 

n=0 \i=Ofe=0 / 



00 00 



(z) = E ^(^ + l)«n+l -2""^ = E + + 2)«n+2 5;". 

n=l n=0 

Consider the particular case 

a{z) = mz + a, b{z) = pz + b, c{z) = qz + c, m,p, q, a,b,c E JR, 

that is the coefficients a{z), b{z),c{z) of the equation (7) are affine functions 
in z. Substituting into (7) gives us 

oo oo 

{mz + a)^{n + 2){n + l)an+2 z^ + {qz + c) E «n+i(n + 1) z^- 

n=0 n=0 

oo / n i \ 

-ipz+b) E E E '^k+i "i-fc+i an-i+i(fc + l)(i - fc + l)(n - i + 1) = 0. 

n=0 \i=Ofe=0 / 

This identity can be written 

oo oo 

m E ''^{f^ + l)an+i -2" + 2aQ;2 + a E (^^ + + l)an+2 -z"— 

n=l n=l 

oo /n— 1 j \ 

~^ E (EE "^=+1 ^i-k+i Oin-i{k + l)(z -k + l)(n - z) J - 6q;?- 

n=l \i=0 fe=0 / 
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oo / n i \ 

~^ X! I X! X] "'fc+i 0!j-fe+i an-i+i{k + \){i-k + \){n-i + \) \ 2;"+ 

n=l \i=Oik=0 / 
00 00 

+q ^ nan + cai + c ^ otn+i in + 1)2;" = 0, 

n=l n=l 

or, equivalent, 

00 

{2aa2 + cai — baf) + ^ [mn(n + l)an+i + a(n + 2)(ra + l)a„+2— 

n=l 

n— 1 « 

XIX! "'=+1 o^n-iik + l){i-k + l){n-i)- 

1=0 k=0 

n i 

-b^Yl "k+i "i-fe+i an-i+i{k + - k + l){n - i + 1)+ 

i=o k=0 

+qnan + c{n + l)a„+i]z" = 0. 
By identifying the coefficients of the powers of z with 0, we find the condition 

2aa2 + cai — baf = 

and the recurrence 

mn{n + l)an+i + a(n + 2)(n + l)a„+2- 

n— 1 i 

-P X X "'=+1 Oin-i{k + - k + l){n - i)- 

i=0 fe=0 

n i 

-6^ X "fe+i '^i-k+i oin-i+i{k + l)(i - A; + l)(n - z + 1)+ 

i=0 A;=0 

+gnQ;„ + c(n + 1)q;„+i = 0, n > 1. (24) 

By the initial conditions 4>{Q) = uq, </)'(0) = ai and 2aa2 + cai — baf = 0, 
this recurrence gives us all the coefficients of the power series (23), but the 
difficult part is just solving the recurrence for the unknown q;„. 

Theorem In the foregoing hypotheses, the multitime series soliton solu- 
tion of the multitime Rayleigh PDE is 

00 

n(a;,t) = X«n(x-AX)", (25) 

n=0 
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with ao, ai fixed, 2aa2 + cai — baf = and an, n > 2 given by the 
recurrence (24) 

mn{n + 1)q:„+i + a(n + 2)(n + l)a„+2— 

n— 1 i 

-P^Yl "fe+i (^i-k+i Oin-i{k + -k + l)(n - i)- 

1=0 fe=0 

n i 

-6^ ^ "fe+i "i-fc+i a„_i+i(fe + l){i-k + l)(n - z + 1)+ 

i=0 A;=0 

+qnan + c{n + l)an+i = Q, n>\. 



6 Families of multitime Rayleigh solitons of 
Van der Pol type 

6.1 Case of coefficients depending on z 

As we have seen in a previous section, if we can choose the constant vector 
Aq. by some conditions which relate the elements h, C, D, A, then we get the 
equation 

a{z)(l)"{z) - d{z)(l)^{z)(l)'{z) +c{z)(p'{z) = 0. 
We can write this equation in the form 

In order to find some solutions of this equation, we make a particular choice 
for Aq,, by a new condition: the coefficients a{z) ^= h°'^{t,x,r],^)XaXi3 — 1, 
diz) = D^{t, x,77,^)A-y and c{z) = C'^{t,x,r],^)X-y will be related by the 
equality 

'a{z)\' c{z) 



d{z) J d{z) ' 
or equivalent 

a'{z)d{z) - a{z)d'{z) = d{z)c{z). (26) 
With such a selection of Aq, the equation becomes 
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that is 

This Bernouhi ODE has separable variables. By integration, we get 

For simplicity, we can take k = and we find a particular family of solutions, 
defined by the relation 



(^) 



If we keep k variable and non-zero, then making the substitution 3k by kf, 
the equality (27) becomes 

Sd(f) f d{z 



f 3# fd{z) , ^TR* 



The integral from the left hand is 

3# 1 /• # 1 /• + 2A;i 



/ci)((/)2 + #1 + A:f) klJ (p-ki kjj + + 

1 , 1^ , I 1 /■ 2(^ + A:i 3 /• # 

= — m |(p — Kil — 



A;? " 2klJ 02 + #i + fc2 ^ 2fciy / 3A:? 



= In |(^ - A;i I - ^ ln((^^ + + ) - "^arctg ^ ^ j + ci , ci G IR. 

Therefore, the general solution of the equation (27) is expressed implicitly 
by the equality 

' m \£zi^ - ^arctg f ^) = / ^ d., e R*. 

Summarizing, wc can formulate the next result: 

Theorem // we can take the constant vector Xa so as the elements 
h°'f^{t,x,r),^), B'f{t,x,r],^), C^{t,x,r],^) and Xa, a, ^,7 = l,...,m, to be 
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related by the conditions 

X, 77, e)AaA^ - 1 = a(x - A„t") 7^ 0, 
BT(t,x,r?,0A-y = 6(x-A«t"), 
C'^{t, X, rj, C)A-y = c{x — Xat°'), 
a'{z)d{z) - a{z)b'{z) = d{z)c{z), 

then the function 

u{x, t) = (f){x — Xat"') 

represents a multitime soliton-solution for the multitime PDE Van der Pol 
(3'), for every (j) defined implicitly by one of the equalities 

^\^) = J_f^^dz (28) 

or 



3 J a{z) 



^1 ^^■'{z) + ^{z)ki + kl ^1 V feiV3 J J a{z) 



(29) 



6.2 Case of constant coefficients 

If we fix the elements h,C,D,X by the conditions 

h'^/^{t, X, rj, C)AqA^ — 1 = constant = a 7^ 0, 
< D^{t, X, rj, 4)A^ = constant = d, 
C'^{t,x,r],^)X-y = constant = c, 

then the ODE (7') takes the form 

acj)" -d(l)'^(f)' + c(j)' = 0. (30) 

Integrating, we find a first order ODE, with separable variables 

acf)' -^(j)^ + c4) = k, kelR. 

In order to find some solutions of this equation, we take k = 0. By this 
choice, the equation becomes a Bernoulli ODE. We make a change of the 
unknown function. 
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and the Bernoulli ODE becomes a linear ODE, 

,, 2c , 2d 

V = 

a 6a 



with the solutions 



— dz 
iP^z) = eJ a 



K 



f 2d 
J 3^ 



2d - I —dz 
e J a dz 



, K eM, 



that is 



Since 



2c 

tp{z) = Kea + K eJR. 

oC 



ip 2 , we obtain 

4>{Z): 



1 



2c 

Kea + — 
3c 



(31) 



(32) 



Thus we obtain a family of multitime soliton-solutions of the multitime Van 

der Pol PDE (3') and we can formulate the next theorem: 

Theorem // we take Aq so as to fix h, D,C, X by the conditions 

h'^^{t, X, 7], ^)A„A^ — 1 = constant = a ^ 0, 
D'^{t, X, T], ^)A^ = constant = d, 
C^it, x,r],^)X^ = constant = c, 

then we get a family of multitime soliton-solutions of the multitime Rayleigh 
wave equation of Van der Pol type (3'), 



u{x, t) 



\ 



^ 3(C7A^) 



(33) 
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